Cold collision shift cancelation and inelastic scattering in a Yb optical lattice clock 



A. D. Ludlow, N. D. Lemke, J. A. Sherman, and C. W. Gates 
National Institute of Standards and Technology, 325 Broadway, Boulder, CO 80305, US-^^ 

G. Quemener, J. von Stecher, and A. M. Rey 
JILA, National Institute of Standards and Technology and University of Colorado, Boulder, CO 80309-0440 

Recently, p-wave cold collisions were shown to dominate the density-dependent shift of the clock 
transition frequency in a ^^^Yb optical lattice clock. Here we demonstrate that by operating such 
a system at the proper excitation fraction, the cold collision shift is canceled below the 5 x 10~^* 
fractional frequency level. We report inelastic two-body loss rates for ^Pq -^Pq and -^Po scattering. 
We also measure interaction shifts in an unpolarized atomic sample. Collision measurements for this 
spin-1/2 ^'^^Yb system are relevant for high performance optical clocks as well as strongly-interacting 
systems for quantum information and quantum simulation applications. 



Large ensembles of ultracold atoms offer atomic clocks 
a measurement of the atomic state with high signal-to- 
noise ratio. For clocks utilizing optical transitions, this 
has the potential to yield time and frequency measure- 
ments with new levels of precision and speed (e.g. [Tl[5]). 
However, large ensembles of cold atoms can lead to high 
number density and thus significant interatomic interac- 
tions. These interactions can perturb the clock transi- 
tion frequency, compromising the accuracy of the atomic 
standard. For example, in cesium fountain primary stan- 
dards, cold collision shifts can become significant [3111], 
influencing clock operation (e.g. [SHT]). 

Gptical lattice clocks, which probe the ultra-narrow 
^Sq -^Pq transition in two- valence-electron atoms held in 
an optical potential, are also susceptible to cold colli- 
sions. Non-negligible collision effects have been observed 
in lattice clocks using nuclear-spin-polarized, fermionic 
samples of ^'^Sr 015] and ^'^^Yb [ID]. Measurement and 
control of these collisions therefore play a key role in the 
continued development of these standards. At the same 
time, the control of these interactions are an integral part 
of proposals for quantum information 112j and quan- 
tum simulation of solid-state-analog Hamiltonians [Tbi - 
lisi . For ®^Sr, the primary interaction giving the cold 
collision shift was identified as an s-wave interaction be- 
tween atoms in non-identical superpositions of the clock 
states [HI HH HZ] . It was recently observed that strong 
interactions can isolate and suppress the shift [16, 18 . 

For ^''^Yb, it was shown that a p-wave interaction be- 
tween atoms in the ^5*0 and ■^Pq electronic states was the 
dominant mechanism responsible for the cold collision 
shift [19]. Unlike the s-wave case, such an interaction is 
less sensitive to small particle distinguishability between 
the ultracold fermions. Gonsequently the observed shift 
exhibits roughly a linear dependence on the '^Pq excita- 
tion fraction (Fig. 1(a)). Moreover, in the weakly inter- 
acting regime of the 1-D lattice, the cold collision shift 
crossed zero at a mean excitation fraction close to 0.5. 
In this work, we exploit this zero-crossing to demonstrate 
cancelation of the ^^^Yb cold collision shift in a 1-D opti- 



cal lattice below the 5 x 10~^^ fractional frequency level. 
The cancelation is enabled by the fact that near 50% ex- 
citation, population of the two electronic states of one 
atom induce equal collision shifts on the opposite states 
of another atom, leaving zero net shift for the clock tran- 
sition 20 . The collision shift cancelation can thus be 
likened to the Stark shift cancelation which optical lattice 
clocks exploit by operating at the "magic" wavelength 
pn [35] . Our uncertainty in this cancelation reaches be- 
low the smallest total uncertainty levels reported to date 
for any type of atomic clock, making it a powerful and 
pragmatic technique for mitigating the collisional shift in 
a lattice clock. 

The ^^^Yb (/ — 1/2) optical lattice clock is described 
in detail elsewhere [TOl [19]. After dual laser-cooling 
stages, atoms are loaded into a 1-D optical lattice using 
a laser at Amagic — 759 nm. Atomic population, initially 
split between the two ground magnetic substates, is opti- 
cally pumped to a single spin state using cr-polarized light 
resonant with the % (/ = 1/2) - ^Pi (/ 3/2) transition 
in the presence of a bias magnetic field (5 G). We esti- 
mate the average atomic density can reach po ~ 3 x 10^^ 
atoms/cm^, with an atomic temperature of ^ 10 fiK. 
The atoms are then spectroscopically probed on the ^5*0 - 
■^Po transition using Ramsey spectroscopy. The Ramsey 
pulse times used here are 5-10 ms, while the dark time 
is 80-150 ms. A cavity-stabilized probe laser is actively 
stabilized to the center Ramsey fringe using an acousto- 
optic modulator. The cold collision shift is the measured 
frequency difference between interleaved atomic samples 
of high and low density. 

Measurement of the collision shift in a 1-D lattice is 
shown as filled points in Fig. 1(a), as a function of ex- 
citation fraction (during the Ramsey dark time). The 
theoretical calculation of the shift using a p-wave model 
is also shown (solid line). In the model, the domi- 
nant interaction, Veg, is between ^5*0 and ^Pq atoms, and 
a weaker interaction between two excited state atoms, 
Vee = 0.1 X Veg, is also included [19]. The shift has a 
zero value near 50% excitation, and this zero-crossing is 
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FIG. 1: (color online) (a) Cold collision shift as a function 
of excitation fraction, for a polarized sample in a 1-D opti- 
cal lattice oriented vertically (filled triangles) or horizontally 
(filled circles). Open squares are with an unpolarized atomic 
sample, (b) Cold collision shift cancelation: the measurement 
uncertainty for this cancelation is shown, given by the total 
deviation (similar to two-sample Allan deviation [23]). (c) 
Measurement of the residual cold collision shift when operat- 
ing at 51% excitation and an atomic density of pi. The red 
dashed line indicates the weighted mean of the seven mea- 
surements, +2.5 mHz, while the solid lines gives the one-rr 
error bars of ±2.4 mHz. 



the focus of our attention here. To determine the precise 
excitation corresponding to zero shift, we made real-time 
measurements of the excitation fraction by turning off 
the second Ramsey pulse and measuring the ground and 
excited atomic populations after the dark time. These 
measurements were interspersed between cycles of usual 
Ramsey spectroscopy. We then slightly adjusted the 
probe laser power to keep the measured excitation frac- 
tion constant during each shift measurement. 

The precision of the zero-crossing measurement is given 
by the clock instability during interleaved measurements 
of high and low atomic density. Our measurement sta- 
bility benefits from recent improvements to the cavity- 
stabilized laser used to probe the clock transition [1]. 
Fig. 1(b) shows the precision of one such measurement, 
where we observed a collision shift of -1.4 (3) mHz after 
14500 s of averaging. Fig. 1(c) shows the result of seven 
similar measurements, taken sequentially over the course 
of several weeks and with different measurement dura- 
tions. For all but one measurement, the mean excitation 
fraction was controlled to 51±0.3±1.3 %, where the first 
uncertainty is the fluctuation in the measured excitation 
fraction and the second is the systematic uncertainty in 



the absolute excitation value. (For data point number 
two, the excitation fraction was set to 1% higher; we 
thus applied a —5 mHz correction to the measured col- 
lision shift, as determined from the slope of the curve 
in Fig. 1(a).) The weighted mean of the seven measure- 
ments (red Hues) is 2.5 (2.4) mHz (reduced = 1.04). 
This corresponds to a fractional shift of 4.8 (4.6) xlO~^^ 
of the transition frequency, and demonstrates the small- 
est measurement of a collision shift in a lattice clock. 

In order to routinely implement this collision shift can- 
celation, we consider the robustness of this technique 
to relevant experimental conditions. For an operational 
density of pi ~ 3 x 10^" atoms/cm'^, a 1% change in the 
excitation fraction leads to a change in collision shift of 

1 X 10^^''. As described above, it is straightforward 
to keep the excitation fixed at or below the 1% level. 
However, a complication arises from time-dependent ex- 
citation due to trap loss. Particularly, we have observed 
inelastic two-body losses involving both '^Pq-^Pq and ^Sq - 
•^Po interactions. With both ^Sq {g) and ^Pq (e) popula- 
tions present, the number density rate equations are: 

fle{t) = -jSeeTleitY - T eTleit) - I3egng{t)ne{t) 

For a single population (ug or rie) Peg loss can be ignored. 
Fig. 2(a) shows ^5*0 trap loss (black triangles, for a pure 
sample of ^Sq atoms) , with a one body loss fit yielding a 
trap lifetime of l/Fg = 480 (20) ms. Also shown is ^Pq 
trap loss (blue circles, for a pure sample of ^Pq atoms). 
This population experiences notably stronger decay at 
high densities, and good fit requires both one-body (Fg) 
and two-body (/3ee) losses. After integrating the losses 
over the spatial extent of a single lattice site, as well as 
across the distribution of occupied sites, we find 1/Fe = 
520 (28) ms and /3ee = 5 x 10^^^ cm^/s, somewhat larger 
than the "^Pq-^Po decay measured in ^^Sr [211 dS]. Because 
the absolute atomic density is difficult to calibrate, the 
uncertainty in Pee is estimated as 60%. 

Inelastic two-body loss for "^Pq atoms can be consid- 
ered with a single atom-atom channel, time-independent 
quantum formalism similar to |261 12 7j . Lacking an accu- 
rate Yb2 potential, the short range physics is described 
here by a boundary condition at an interatomic sep- 
aration Pq = 20 ag (ao is the Bohr radius), repre- 
sented by two parameters: (i) an accumulated phase 
shift, < 5 < 7r, due to the unknown short-range po- 
tential from R — to R — Rq, and (ii) a loss probability, 
< Pis < 1 at P = Pq- The long range physics is ac- 
curately described by a van der Waals interaction, with 
Cg = 3886 a.u. given by [5S]. The Schrodinger equa- 
tion is numerically solved from Rq to P — >■ oo, and the 
cross section is computed as a function of the collision 
energy. The thermalized rate coefficients are found by 
averaging the cross sections with a Maxwell-Boltzmann 
distribution of collision energies for a given temperature. 
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FIG. 2: (color online) (a) ^So trap loss (black triangles) with 
an exponential fit (dashed curve). "^Pq trap loss (blue circles) 
with a one- and two-body loss fit (solid curve). The inset 
shows loss for a pure sample (black triangles) and a 
mixed sample of ^S'o and ^Po atoms (blue circles) . The former 
is fit with a simple exponential, the latter with ^Pq -^Pq and 
^So -^Po two body losses, (b) Two body loss rates as a function 
of temperature, calculated with pig = 1 (dashed curves) and 
with Pis = 0.8 and 5 = 0.51 tt (solid curves). Black is for 
a nuclear-spin-polarized sample of ^Po atoms (/3ee), red an 
unpolarized sample of ^Po atoms ([/See + /3ee]/2). Circles give 
experimental measurements, black for /See and blue for /Jeg. 



With full loss probability (pig = 1), /3ee is shown in Fig. 
2(b) (black dashed curve) as a function of temperature. 
Here, the rates are universal, independent of 6 [S^ WJ\ . 
The interaction involves two identical fermions, and thus 
quantum statistics dictates that the interactions must be 
odd partial waves, notably p-wave at these ultracold tem- 
peratures. As two excited atoms approach each other, 
those which successfully tunnel through the p-wave bar- 
rier are assumed to inelastically scatter with unit prob- 
ability, a reasonable assumption due to the large num- 
ber of exit collision channels available to a pair of atoms 
with high internal energy. For an unpolarized sample 
of ^Pq atoms (here with equal to/ = ±1/2 populations), 
interactions between distinguishable atoms also include 
an s-wave term, and the loss rate for distinguishable ^Pq 
atoms is labeled /3ee. The total loss rate for the unpo- 
larized sample, given by the average of /3ee and /3eej is 
shown in Fig. 2(b) as a red dashed curve (pis = 1). Ex- 
perimentally, we observed identical loss rates within 10% 
for a polarized and an unpolarized sample at 10 //K (Fig. 
2(b), black circle). Since full loss predicts unequal loss 
rates for polarized and unpolarized samples at 10 //K, ^Pq 
atoms may not be lost with full unit probability at short 
range. Instead, using short range parameters of pis = 0.8 



and (5 = 0.51 tt, the loss rates are shown as solid curves 
in Fig. 2(b) for polarized (black) and unpolarized (red) 
cases. Here, we found better agreement with the exper- 
imental data, suggesting a deviation from the universal 
regime. 

For a mixed population of and ^Pq atoms, we ob- 
served additional loss through /3eg. Fig. 2(a) inset high- 
lights this by comparing ^Sq trap loss for two different 
atomic samples: with (blue circles) and without (black 
triangles) the presence of ^Pq atoms. The additional loss 
is particularly notable at short times where both ^S'o and 
•^Po populations are large and is consistent with a ^5*0 -^Pq 
two-body loss rate at the level of /3eg = 3 x 10~^^ cm'^/s. 
The ^Pq population is excited coherently from the ^Sq 
state, but due to excitation inhomogeneity, the lossy col- 
lisions are between partially distinguishable atoms [TO] . 
The magnitude of the inelastic loss is noteworthy because 
for molecular states correlating to a ^Sq and '^Pq atom 
pair, only the ground state (correlating to the ^Sq- 
^5*0 state) lies at lower energy, and long-range coupling 
to this state is spin-forbidden. 

In general, all of these loss processes lead to a time- 
dependent excitation fraction during the Ramsey dark 
time. This, in turn, affects the balance between -^^o and 
•^Po collisionally induced energy shifts. For simplicity, we 
can easily operate at a lower atom number density, reduc- 
ing not only the collisionally-induced shifts, but also the 
two-body inelastic losses. At the same time, the number 
of quantum absorbers remains well in excess of 1000, in 
order to accommodate a high signal-to-noise ratio. At 
a density of pi, the excitation fraction over a dark time 
of r = 150 ms changes by only several percent. We 
therefore take the time-averaged value to indicate the 
excitation fraction where the shift is canceled. 

Any degree of imperfect polarization of the nuclear 
spin state introduces a host of other possible atomic inter- 
actions. In p-wave, the Vg" interaction, which shifts the 
singlet state, becomes allowed [in]. Furthermore, distin- 
guishability between different nuclear spin states allows 
s-wave interactions C/gg, J7ee, and . In -^''-^Yb, both the 
s- and p-wave \gg) interaction terms are small |29j . How- 
ever, the remaining interactions can alter the observed 
collision shift from the spin-polarized case. For this rea- 
son, high polarization purity is important for optical lat- 
tice clocks. Here, we benefit from the simple structure 
(/ — 1/2) of -'^^^Yb and can readily optically pump to a 
single spin ground state with 99% purity. The polariza- 
tion purity is directly measured by observing the absence 
of a clock excitation spectrum for the unpopulated spin 
state. 

To quantify the effect of imperfect polarization, we 
measure the collision shift as a function of excitation for 
unpolarized atoms (open squares in Fig. 1(a), equal mix- 
ture of both spin states). Since, during spectroscopy, 
a weak bias field (1 G) lifts the degeneracy of the tt- 
transitions from the two spin states, here the excitation 
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FIG. 3: (color online) Cold collision shift as a function of 
excitation for a polarized sample in a 2-D lattice, for different 
Ramsey times. p2 is a typical 2-D lattice density, where 25% 
of the atoms are doubly occupied in a lattice site with density 
of 4 X lO^'^ atoms/cm'^. Circles correspond to experimental 
measurements, and solid curves are theoretical calculations 
|19| . In both cases, blue is a Ramsey time of T = 10 ms, red 
is 40 ms, black is 80 ms, green is 160 ms, and orange is 210 
ms. 

fraction is defined as that for the particular spin state 
being resonantly excited, not accounting for the other 
'spectator' spin state. In general, the measured shifts are 
smaller than for the spin-polarized case, implying that 
competing interactions have the opposite sign as those in 
the polarized case. Notably, the measured zero-crossing 
in the shift occurs at a lower excitation fraction (around 
40%), leaving a net positive shift at 51%, where the shift 
is zero for the spin polarized case. Based upon these 
measurements, we determine a 1% polarization impurity 
will not affect the shift zero-crossing at or above the level 
measured in Fig. 1(c). 

We have shown that the cold collision shift can be can- 
celed at the 5 x I0~^^ level in a 1-D optical lattice clock of 
^^^Yb. Other implementations may be suitable for reduc- 
ing the clock shift at or below this level. In a 2-D optical 
lattice clock, strong interactions can lead to a decay of 
the collision shift [19j|30]. In particular, longer Ramsey 
dark times lead to smaller shifts and can reduce the shift 
dependence on excitation fraction (Fig. 3), making it at- 
tractive for shift reduction. Care must be taken since, as 
we have observed both experimentally and theoretically, 
the interactions can also reduce Ramsey fringe contrast. 
Alternatively, the 2-D lattice system also exhibits a zero 
crossing in the shift versus excitation. However, strong 
interactions can move the zero-crossing excitation with 
a nonlinear dependence on interaction strength (Fig. 3), 
and thus the weak interactions of the 1-D lattice may be 
easier to control. 

The 3-D optical lattice continues to be an interesting 
choice [2Tj, where it is straightforward to achieve high 
atom number with an average of <C 1 atom per lattice 
site. The small fraction of lattice sites with double oc- 
cupancy will exhibit very strong interactions. Doubly- 
occupied sites could also be eliminated using photoas- 



sociation losses [3T] or directly via the two-body losses 
observed here. Furthermore, the 3-D optical lattice sys- 
tem may exhibit kinematic suppression of the atomic in- 
teractions which yield the shifts ■ Indeed, these rich 
atomic systems will undoubtedly continue to offer many 
interesting phenomena in 1-, 2-, and 3-D confinement. 
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